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ABSTRACT 


I A characterization  of  the  best  L^ -approximation  to  a continuous  function 

I by  classes  of  fixed-knot  polynomial  splines  which  satisfy  generalized  convexity 

constraints  is  presented  emd  uniqueness  is  shown.  Included  is  the  possibility  of 
specifying  the  positivity,  monotonicity,  or  convexity  of  the  class.  The  proof  of 
uniqueness  uses  recently  developed  results  for  Hermite-Birkhoff  interpolation  by 
^ splines. 
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SIGNIFICANCE  AND  EXPUU^ATION 


In  practical  applications  of  approxisation  theory,  it  often  occurs  that  the 
function  being  approximated  is  known  to  have  sosm  additional  properties  such  as 
positivity,  mono tonicity,  and/or  convexity.  A best  approximation  may  not  preserve 
any  of  these,  even  though  such  properties  may  be  very  impoxtant  for  the  application. 

if 

\ In  this  case  one  should  search  only  astong  approximations  which  satisfy  the  desired 

£ 

i side  conditions  of  positivity,  sonotonicity,  or  convexity.  These  are  special 

j cases  of  generalized  convexity  constraints. 

^ j This  paper  begins  the  study  of  best  approximation  in  the  mean  or  L^^-norm  by 

fixed~knot  polynosiial  spline  functior^  which  satisfy  generalized  convexity  con- 
straints. As  one  particular  example,  it  is  sho%m  that  the  best  L^-approximation 
to  any  continuous  function  by  nonnegative,  nonincreasing  cubic  splines  with  fixed 
knots  exists  and  is  unique. 
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The  responsibility  for  the  wordinq  and  views  expressed  in  this  descriptive  summary 
lies  with  MRC,  and  not  with  the  author  of  this  report. 
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BEST  MEAN  APPROXIMATION  BY  SPLINES  SATISFYING 
GENERALIZED  CONVEXITY  CONSTRAINTS 

Dennis  D.  Pence 


Introduction 

The  concept  of  monotone  approximation  by  polynoetials  was  introduced  by  O.  Shisha  1121 . 
It  has  been  studied  by  many  authors  including  R.  A.  Lorentz  |4),  who  demonstrated  unique- 
ness for  best  approximation  by  scnotone  polynomials  in  the  uniform  and  L^-norm.  J.  A. 
Roulier  Md  G.  D.  Taylor  t9]  have  generalized  the  concept  to  include  sore  arbitrary  re- 
strictions on  the  ranges  of  derivatives.  If  the  restrictions  are  all  of  nonnegativity  or 
nonpositivity,  they  are  called  generalized  convexity  constraints. 

Classes  of  splines  satisfying  generalized  convexity  constraints  were  introduced  and 
studied  in  the  author's  thesis  (6).  Best  uniform  approximations  were  characterized  and 
partial  uniqueness  was  established.  This  paper  continues  the  study  of  such  splines  by  con- 
sidering the  L^-norm.  A characterization  for  best  L^-approximation  by  these  classes  of 
splines  is  given  and  uniqueness  is  demonstrated.  Previous  papers  considering  best  L^- 
approximation  by  splines  include  the  work  of  M.  P.  Carroll  and  D.  Braess  (2) , R.  V.  Galkin 
[3) , and  A.  Pinkus  |8| . 

we  use  results  f>r  Hermite-Birkhof f interpolation  by  splines  which  were  developed  by 
the  author  [6],  (71.  For  cos^leteness  the  required  results  are  given  in  the  following 
section. 

1.  Hermite-Birkhof f Interpolation  by  Splines 

Suppose  -<»<aE£  ,=  b<‘"  and  integers  R,  with 

" 0 1 q q+1  v 

0<R  <m,  V”l,...,q,  are  given.  Let  S"  ■ 5”*({C  {r  denote  the  space  of 

V—  PPVj^Vj^ 

polynomial  spline  functions  of  order  m with  fixed  knots  each  with  multiplicity 

Q _ 

R , respectively,  where  p = ^ R . Thus  g e S is  piecewise  a polynomial  of  degree 

V V P 


at  most  m - 1 with  g 


(j) 


u-1 


discontinuous  only  at  a knot  ^ 


where  j > m - R . We 
— V 


adopt  the  convention  that  all  elements  of  and  all  derivatives  of  elements  of  ^ are 
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defined  everywhere  by  assuninq  continuity  from  the  right.  Notice  that  dims  • m * p. 

P 

We  review  needed  facts  about  Hermte-Birkhoff  interpolation  (HBI)  by  polynoaUal 
splines. 

Let  Interpolation  points 


be  given.  A matrix 


X - {a  < Xj^  < Xj  <•••<  Xj^  < b} 


E - i ■ » j - 0,  l,...,m  - 1 


is  called  a spline  incidence  matrix  for  X and  5^  provided  e^^^  • 0/  1 1>  or  2 and 

e,  . " - 1 or  2 only  if  x.  “ £ for  some  v and  j > m - R . The  HBI  problem  defined 
ij  i 'v  - V 

by  (E,  X,  5^)  is: 

Given  any  values  ' ®^j  • 1 or  2}  and 

(Yij”  ! e^^  « -1  or  2),  find  9 c with 

(1.3)  g*^' (Xj^)  - whenever  e^^  - 1 or  2 , 

(1.4)  g*^’(x^-)  - Y^j”  vdienever  e^^  - -1  or  2 . 

As  in  (6],  (7),  when  we  display  such  a matrix  E,  we  indicate  the  relationship  be- 
_yfeen  the  interpolation  points  X and  the  )cnots  of  the  spline  space  ^ by  drawing  the 
following  lines: 

(i)  If  x^  < C < a solid  line  between  the  1-th  and  (i  + l)-th 

rows  extending  from  the  (m  - R )-th  column  to  the  (m  - l)-th  column.  If 

more  than  one  )cnot  lies  between  x^  and  *14.1'  then  draw  several  lines. 

(li)  If  x^  “ 5 , we  enclose  in  a box  the  entries  in  the  i-th  row  from  the 

(m  - R )-th  column  to  the  (m  - l)-th  column. 

V 

Thus  an  entry  of  E may  be  -1  or  2 only  if  it  is  boxed. 
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In  particular  when  (E,  X.  S^)  la  full,  the  (LPC)  iaply  that 


(1.11) 


LEMHA  1.1.  If 


|E(0  I ».q*l)||  » [ * , I ■ I, 


vl 


B(0  I t.a) I 


♦ p(0  I t,s)  for  aoae  0<I<a£q^l  or 


0<t<8<q>l,  then  (B.  X,  5^  can  be  split  vertically  into  or  three  HBl  pro- 
blem. each  defined  on  a spline  space  of  order  still  ■ but  with  fewer  lutots  than 
Ihe  "central*  one  of  the  decoes>osed  probleaw  has  incidence  sistrix  E(0  i t,s) . 

I for  all  1 ■ 0.1,..., 


LEMMA  1.2.  If  for  s 


I and  V,  and  e^^ 


- » - I, 

V 


or  if  - n l2£.  soMs  v,  then  (B.  X,  5^  can  also  be  split  vertically  Into  two  HBI 
problesM  considering  fewer  knots.  E(0  i 0,w)  and  B(0  t v.q'^l)  will  t>>e  incidence 
■atrices  for  these  two  ssaller  probleas. 

We  further  note  that  the  above  decoapositions  preserve  the  (LPC)  and  t)»at  if  the 

original  problem  has  a full  matrix,  then  so  do  all  of  t)>e  smaller  problems,  (^asi-polsed- 

ness  of  (B,  X,  S*)  is  equivalent  to  quasi-polsedness  of  all  of  thm  split  problems. 

P 

SisUlar  decompositions  have  been  noted  by  several  auttiors,  see  (5),  for  example.  The  com- 
plete details  are  tedious  but  not  hard  and  can  be  found  in  [6] . 

Let  (E,  X,  S^)  indicate  a given  HBI  problem.  If  x^  / '**  tV»t 

we  have  a regular  sequence  beginning  with  e^^  of  order  p when 

i,jtp-l-'  “i.J-l-°  •it''” 

Also  if  t » then  we  say  that  we  have  a regular  sequence  beginning  with  e^^^  of 

order  p when  e^^  - ®i,J+p-l  i "ith  j ♦ K 1 " ' *1. j-i  ’ » 

e . « 0 if  either  is  defined.  Further  a regular  sequence  e ,...,e  is  called 

i,j+p  iJ  ‘■ 

strongly  regular  if  is  defined,  zero,  and,  in  tlie  case  where 

x^«5,j+p<m-R.  A sequence  is  even  if  it  has  even  order  and  odd  otherwise. 

We  say  that  a regular  sequence  e^y  ....e^^  J+P“l  supported  provided  there  exist 
integers  ij^,  j^,  i^,  Jj  with  i^  < i < i^,  e^  ^ • 1 or  2 , 
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< mnlj  , {n-R^  IX.  < 4^  « x^)|  , 


h * ‘ “l  * ^ * “ij»'  . .nd 


r 1,  if  X / {£  >7  . 

J 

(1.14)  o,  • / 1,  if  X,  ■ C xixl  j,  < ■-«  » 

ij.Jj  ij  V 2 V 

-1  or  2,  if  X,  • C and  j.  > m-R 
\ ij  V 2 — V 

The  problem  (E.  X,  5^)  is  called  weaXly  conservative  (C)  if  every  supported  stronqly 
regular  sequence  is  even. 

THEOREai  1.1.  Suppose  (E,  X,  satisfies  the  (LPC)  and  (C) . Then  it  ^ quasi-poised. 

This  theorem  generalizes  the  sufficiency  theorem  of  Atkinson  and  Sha..!Ba  for  HBI  by 
polynomials  (1).  The  proof  can  be  found  in  |6) , (7). 

We  shall  need  the  following  technical  lemnas.  All  three  lemmas  concern  attempts  to 
add  conditions  of  some  sort  to  a given  HBI  problem. 


LEMMA  1.3.  S 


(E,  X,  5*)  satisfies  the  (LPC)  but  when  some  strongl 


sequence  ^ E extended  to  have  an  additional  one  to  the  right  giving  the  matrix  E 
then  the  (LPC)  are  violated.  There  exists  C / (C  )?•'.’  X so  that  when  C is  added  as 


a ^m^le  knot  to  the  spline  space,  then  (E  , X,  5^.,  ((4  )?.  4;  (R  )?,  D)  satisfies  the  (LPC). 

Vi  VI  ■ — " 


Proof.  Suppose  e. 


e,  . , • 1 is  the  strongly  regular  sequence  of  E and  that 
1,  j"i 


e.  . is  changed  from  a zero  to  a one  to  obtain  E.  Then  anil  ,0^i<s<_q  + l 


X,  £ 14,, 4 1 so  that  for  s"  , 
1 £ S p 


(1-15)  ||E(n  : £pS)l|  ■ ||E(n  ! ••  1 

» m - n + p(n  I t»s)  + 1 . 

without  loss  of  generality  assume  that  (1.15)  cannot  happen  first  for  any  q > n sod 
secondly  with  n for  any  I and  s with  I4j'4gl  ^ 
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Let 


(1.16)  e - min  ( • (x^^j-x^) . ( * *|H  • 

Vie  choose  ( c (x^  - c.  x^  > e)  \ (x^)  in  such  a way  that  when  ( is  added  to  the  knot 
set.  i.e.  properly  ordered,  C ■ Cj  . Rj  ■ 1#  knd  t < i < a « 1, 

we  have 

(1.17)  ||E(n  I t.t)||  < ||E(n  t l.s«l)||  and 

(1.18)  ||E(n  : t,s>l)||  < ||E(n  > t.s«l)||  . 

with  respect  to  S^,  ((C)?^^  » {R')?*^) 

P^i  V 1 Vi 

It  is  easily  seen  t)tat  this  can  be  done.  The  proof  is  coapleted  by  checkinq  the 

various  ways  i might  violate  tne  (I.PC)  with  respect  to  tt»e  new  spline  apace.  If  this 

happens,  then  (E,  X,  S^)  must  violate  tive  (LK) , contrary  to  hypothesis. 

P 

A condition  corresponding  to  e,  “ 1 is  called  a Lagrange  condition  and  we  say 

t $0  ' 

that  we  are  adding  a Lagrange  condition  at  t to  an  HBI  problem  (E,  X,  5^)  if  a sero 
in  the  j ~ 0 column  is  changed  to  a one,  possibly  by  adding  a new  row  to  E if  t / X. 
LEMMA  1.4.  Suppose  (E,  X,  £^)  satisfies  the  (LPC)  b^  ^ not  full.  Then  there  exists 
a point  t / X U *2.  that  a Lagrange  condition  can  )ae  added  to  (E,  X,  S^)  with- 

out violating  the  (LPC) . 

Proof.  If  ||e(0  : l,s)  ||  - m + p(0  i t,s)  for  some  0<t<s<q*l  or 
0<f<s<q+l,  then  we  can  decompose  according  to  Lemma  1.1  and  consider  one  of  the 
split  problems  which  is  not  full.  Thus  witlvout  loss  of  generality  we  assuse  this  never 
happens.  But  then  any  Lagrange  condition  can  be  added  wit)x>ut  violating  tlie  (LPC). 

LEMMA  1.5.  Assume  that  (E,  X,  ^)  satisfies  the  (LPC)  but  is  not  full.  Without  loss 
of  generality,  assume  that  (b})  C x,  possibly  ^ having  some  rows  with  all 

zero  entries  ^ E.  Then  we  can  "fill*  E in  such  a that  the  (LPC)  remain  valid  ^ 

(i)  chemqinq  some  "boxed"  zeros  to  ones, 

(ii)  changing  some  "boxed"  ones  to  twos,  and/or 

(iii)  changing  some  zeros  to  ones  in  the  last  row  corresponding  to  b. 


Proof.  Inductively  for  1 • 1,  2.  , q we  sake  chanqec  of  type  Cl)  or  (ii>  for  “bojied* 

entries  correspondinq  to  the  interpolation  point  and  knot  so  that  (1.11)  will  be  valid 

for  that  inteqer  after  the  changes  are  wade.  Further  we  wake  changes  one  at  a tiiw  for 
entries  with  j-index  as  large  as  possible  without  violating  the  (LPC) . show  that  this 
is  always  possiile,  suppose  we  have  done  this  for  t ■ 1.2,...,  ^«~1  (If  any)  and  )>ave 

1 

(1.19)  '|e||  - I|e(0  I l.q*l)I|  > I R , f - 1.2... ..*,-1  . 

” v-1  '' 

Suppose  X “ C If 

(1.20)  ||e||  - ||e(0  1 »..q4l)l|  I I \ . 

v-1 

then  there  is  no  need  to  make  any  changes  in  t)ie  i,-th  row.  If  e.  - -1  or  2 for  all 

j - m - R m - 1,  then  (1.20)  holds.  Suppose 

^ • 

(1.21)  e • 0 or  1 where  ■ - R^  J,  , 

(1.22)  e ■ -1  or  2 for  all  j • (»  - R - 1 (if  any)  • 

1,3  1, 

but  e . cannot  be  changed  to  -1  or  2,  respectively,  wit)x>ut  violating  the  (LPC).  By 

i,.3, 

(1.10),  there  exist  integers  f\,  i,  and  s with  I < 1,  i s,  n ^ j,  . R^  * b - h for  all 
? < u < s (if  ^u^y) , and 

(1.23)  llE(n  : i,s)l|  - m - n + p(n  : *,8)  . 

Then 

(1.24)  IIeII  - I|e(0  : f,,q+l)ll 

> ( IIeII  - I|e(0  1 i,q+l)  II  ) 

C ||E(n  s l,,s)  II  , if  t, 

^n\ax(0,  m-R  -nl  » 


♦ ||E(n  : t,5)  II 

\ ♦ max (0,  n - m ♦ R,  ) 
■»  / 
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I SIBIMI  I II  ~ 


r—T 


» I R . 
v»l 

Again  (1.20)  holds.  Thus  we  can  always  accoaplish  the  induction  step. 

If  the  matrix  is  still  not  full  after  all  of  these  changes,  then  we  make  changes  of 
type  (lil)  at  entries  with  j-index  as  large  as  possible  without  violating  the  (LPC) . Me 
argue  in  a similar  manner  that  if  it  is  not  possible  to  change  some  such  entry,  then  it  is 
unnecessary  to  do  so. 

Example  1.1. 


becomes  via  the  procedure  given  for  Lenma  1.5 


The  display  E is  quasipoised  and  E is  poised  by  Theorem  1.1. 

2.  Best  L^-approximation  by  Splines  with  Generalized  Convexity  Constraints 

Let  integers  0_fk|j  <)Cj^  < •••  1 and  e^“±l,v"0,l,...,w  be  givsn. 

Suppose  C[a,bl,  i.e.  < m,  v “ l,...,q  and  ra>  1. 


Define 

(2.1) 


((c  ) 

G - {*  e Sp  : (t)  > 0 


a < t < b , V « 0,1, 


w) 


Recall  we  assume  right  continuity  of  all  spline  derivatives.  Also  for  every  g e G we 
have 


(2.2) 


g (C„-)  L 0 if  ^ m 


-8- 


For  an  integrable  function  h.  let  ||h||^  - /^|h(t)  |dt  . Suppose  f is  in 
C(a,b],  the  space  ot  continuous  functions  defined  on  (a,bl , but  is  not  in  G.  Then  a 
best  L^-approximation  to  f frosi  G is  a spline  g,  t G such  that 


(2.3) 


|g.  - f|lj  “ inf  ||g  - f|lj 
G 


Denote  by  P^Cf)  the  collection  of  all  such  best  approximations.  P(.(^)  • ♦ because  G 
is  closed,  convex,  and  finite-diuensional.  We  have  the  following  characterisation  theorem. 
Assume  f c C[a,b]  and  G is  defined  as  in  (2.1).  Tben  there  exist 

(i)  functions  , r > 1 where  |v>^(t)  | - 1 for  almost  every 

t e la,b)  , i « 1, . . . ,r  , 

(ii)  an  HBI  problem  (E,  X,  5^)  , ||e||  + r£m  + p + l where  e^^  * o only  if 

j ” for  some  v, 

(iii)  positive  scalars  and 

(iv)  scalars  : e^^  - 1 or  2)  , (u.^"  : e^^  - -1  or  2)  where 

sgn  11^^  • - sgn  whenever  j » 


with 


(2.4) 


r D 

^ (t)<-(t)dt  + (x.)  ! e.  . - 1 or  2) 

i.l  1 a 1 ij  i i] 

+ ♦^^’(x^-)  I e^^  - -1  or  2}  - 0 , 


for  all  ill  e , 

p 


such  that  g £ G i£  in  P^(f)  if  and  only  if 


(2.5) 


(2.6) 


(2.7) 


(t)  [g(t)  - f(t))dt  = ||g  - f||  , i - 1 r , 

(j), 


g"  (Xj^)  « 0,  whenever  e^^  » 1 or  2,  euid 


(Xj^-)  ” 0,  whenever  e^^^  » -1  or  2 
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Proof.  This  theorem  is  a specific  case  of  a theorem  of  Rozema  and  Smith  llO,  Theorem  4.1) 


once  we  note  that  it  is  easy  to  find  a polynomial  ^ satisfying 

(2.8)  ♦ (t)  > 0 ) a < t < b;  V - 0,l,....w  . 

LEMMA  2.1.  ^ f c C[a,b]  , qj^  and  g^  are  elements  of  PQ(t) . end  g^^  • ♦ 92*  * 

then  g^^  - vanishes  at  the  zeros  of  g^  - f . 

This  lemna  can  be  found  in  [2]  and  is  a special  case  of  (10,  l^sma  6.1). 

TOEOREM  2.2.  For  every  f < C[a,b),  there  is  £ unique  best  L^-approximation  from  C 
defined  as  in  (2.1). 

Proof.  Let  (E,  X,  ^)  , «n<l  be  as  guaranteed  by 

Theorem  2.1.  Without  loss  of  generality  we  may  assune  (E,  X,  S^)  is  quasl-)poised  since 
any  dependency  in  these  conditions  could  be  used  to  accomplish  (2.4),  (2.6),  and  (2.7) 
with  a smaller  HBI  problem  made  up  of  independent  conditions. 

Suppose  g^^  and  g^  are  both  in  • Then  g^  ■ 'j(9j^  ♦ 92* 

since  it  is  convex.  Applying  (2.5)  to  g^,  we  see  that  alisost  everywhere  that 
gg(t)  * f(t)  we  have 

(2.9)  f^(t)  = sgn(gg(t)  - tit)],  i « l,...,r  . 


Let  t be  the  right  continuous  function  defined  by 


(2.10) 

Define 


♦(t)  “ lim  sgn(g  (t+c)  - f(t+E)] 
e+0 


(2.11)  T = closure  (t  £ (a,b]  : ♦(t)  = 0 or  ♦ changes  sign  at  t)  . 

Thus  9Qtt)  = f(t),  for  every  t € T (an  possibly  at  other  isolated  points  where  q^  - 
doesn't  change  sign) . 

Suppose  g € P_(f).  If  e = 1,  then  J » )t  for  some  v and  g*^Nx.  ) « 0. 

G ij  V 1 

By  the  definition  of  G in  (2.1),  g^^*  does  not  change  sign.  Using  the  zero  counting 

procedure  for  splines  devised  by  SchumaJter  (111  (see  also  (61,  (71),  this  means  that  x^ 
is  in  some  interval  (possibly  degenerate)  where  g*^^  is  identically  zero  and  that  this 


-10- 


interval  is  either  an  even  zero  for 
elude  that  if  a<x^<b,  j<tn-l, 


or  it  contains  one  of  (a)  or  (b).  We  con- 

■ 1,  and  e.  . , is  not  “boxed*,  then 
1 # 


(2.12) 


^(j*i) 


(x^)  • 0 . 


One  by  one,  we  add  enough  of  the  zero  conditions  designated  in  (2.12)  to  the  HBI  pro- 

blem  (E,  X,  Op)  to  assure  that  it  has  no  odd  strongly  regular  sequences  in  rows  for 

which  a < X.  < b and,  if  the  sequence  begins  with  e.  , for  which  x,  / T.  We  add  con- 

ditions  from  (2.12)  to  (E,  X,  •'^)  . if  necessarv,  to  assure  that  it  has  no  even  strongly 

regular  sequences  beginning  with  e if  x.  r T (a,b)  . If  it  is  necessary  at  each  step 

1 « u i 

in  order  to  preserve  the  (LPC) , a simple  knot  is  added  to  the  spline  space  as  described  in 

Lemma  1.3.  Then  we  define  scalars  u.  .,  possibly  zero  and  having  no  particular  sign,  so 

^ . J 

that  (2.4)  remains  valid,  even  for  the  enlarged  spline  space.  The  fact  that  such  constants 
^i  j can  always  be  chosen  follows  from  elementary  linear  algebra  and  Theorem  1.1.  Thus  we 
obtain  (E,  X,  5^)  < and  satisfying  the  (LPC)  of  the  desired  form  with 

r b - Ml 

(2.13)  I ^ (x^)  I - 1 or  2) 


+ = ®ij  - -1  or  2)  - 0 , 

for  all  </  « 

P 


and  for  any 

g e r^(f)  C C ^ 

, ( in 

particular 

for 

gj^  and  g^) 

(2.14) 

g'^’(x.) 

- 0, 

whenever 

"ij  “ 

1 or  2,  and 

(2.1^) 

g‘'>(x,-) 

« 0, 

whenever 

"ij 

« -1  or  2 

Choose  a maximal  subset  tj^,...,t^  from  T \ {x^  i that  when  Lagrange 

conditions  at  t, ,...,t  are  added  to  (E,  X,  5^)  to  obtain  (E  , X , S?) , the  (LPC) 
la  P X X 

are  still  satisfied.  Note  that  (Ej^,  X^,  S^)  satisfies  (C)  so  by  Theorem  1.1  it  is 
quasi-poised. 


Claim  1 ! 


m + p 


If  not,  then  by  Lenne  1.4  there  exists  a point  t.  ^ * Lagrange 

condition  at  t^,  can  be  added  to  get  (E^.  X^,  . t,  ^ T or  else  it  viould  have  been 

added  in  before  using  the  maximality  of  the  chosen  tj,...,t^. 

We  obtain  a full  (?,  X,  5^)  by  applying  Lcmna  1.5  to  (E^,  Xj,  S?),  if  necessary. 
The  only  Lagrange  condition  that  might  possibly  be  added  in  this  step  is  one  at  b.  Note 
that  this  full  HBl  problem  will  still  satisfy  (C)  as  well  as  the  (LPC)  so  by  Theorem  1.1 

m 

it  is  poised.  There  exists  a unique  f with 


(2.17) 


' - ♦(t.)  * 0 , 

(x^)  « 0,  whenever  e^^^  ■ 1 or  2 and  * t^  , 

V “ 0,  tdienever  e^^  » -1  or  2 


First  i|),(t)  - 0 for  all  t e T.  If  t e T \ (tj^, . . . ,tjj) , then  t £ 

where  equality  occurred  in  the  (LPC)  indexed  by  (0,  t,  s)  for  ' hence  in 

(E,  X,  5?)  and  t / [i  ,(  }•  (E,  X,  5^)  can  be  decomposed  according  to  Lemma  1.1  and 

p ‘Is  p 

the  split  problem  involving  t will  have  all  zero  data  values  from  (2.17).  Thus  is 

identically  zero  on  ISj'Sgl- 

If  ||e(0  : l,s)  II  = m + p(0  : l,s)  for  some  0_<I<sf^q  + l where  t,,  / (t^.t^l. 
then  the  problem  decomposes  according  to  Lemma  1.1  again.  Examining  the  part  of  (2.17) 
that  each  split  problem  must  satisfy,  we  conclude  that  i|i,(t)  ” 0 for  all  t e (a,CgJ  if 
5^  < t^  or  \))„(t)  “ 0 for  all  t e if 

If  X,  • E * t*  and  e,.  » 1 for  all  j » 0,1,... ,m  - R * 1,  then  (E,  X,  5?) 
i ’v  * ij  v P 

can  be  decomposed  according  to  Lemna  1.2.  As  above  the  split  problem  not  involving  t,, 

will  be  homogeneous  so  i(i,^  is  identically  zero  either  to  the  right  or  the  left  of  x^^. 


Suppose  i|i,(t)  ” 0 for  all 

t.  / (E,  .E„  )•  Let  (E*,  X.  S^) 
o °0 


t £ IE,  .E  1 for  some  0 ^ < s^  < q + 1,  hence 

*■0  ®0 

denote  (E,  X,  s")  with  the  Lagrange  condition  at  t 
P 
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deleted.  Let  t.  be  any  point  in  (C.  >C  ) \ X.  A Leqrange  condition  at  t.  cannot  be 

0 0 

added  to  (E  , X,  5?)  without  violatlnq  the  (LPC) . If  it  could  be,  it  would  give  a full 
P 

poised  HBI  problem  for  which  the  nontrivial  i|i^  satisfies  all  sero  data  values  which  is 
impossible.  Thus  there  exist  integers  0<^t<s^q'fl  with  i ^ t,  i ' 

and  ||e  (0  : 1,8)  ||  - m t p(0  : t,s) , hence  ||e(0  : t,s)  ||  > m -»  p(0  i t.s) . As  before 

either  i|)^(t)  » 0 for  all  t t [a,C  1 if  C < t*  or  <i,(t)  » 0 for  all  t < U,  ,b] 

®o  “o  ‘o 

if  t,  < e,  . 

*0 

We  conclude  that  is  identically  zero  except  on  some  knot  Interval  (C,  1 

containing  t^  (possibly  (a,b]).  On  this  interval  there  are  only  a finite  number  of 
points  from  T and  only  a finite  number  of  points  where  is  zero.  Any  sequence  be- 
ginning with  e^  Q for  which  * *i  * *i  * strongly  regular. 

Further 

(2.18)  ||e  (0  : l,s)  II  < m + p(0  : f,s)  if  0£f<S£q  + l 

and  (?,.e„)  ^ (C,  ♦ . 

Claim  2:  son  i(i*(t)  - «(t)  , for  a.e.  t e (£.  .C  I 

Suppose  ijij,  changes  sign  at  some  t e (C,  iC  ) \ (T  U {t  })  where  ♦ does  not 
change  sign.  If  t « x.  for  some  i and  e”.  =1,  then  e,  begins  a strongly 

Ux  XfU  XfO 

regular  even  sequence  e , . . . ,e  and  since  t must  be  an  odd  zero  for 

X # u X f y*x  0 

I)!*,  = 0.  Define  (E^,  to  be  (E  , X,  S?)  with  e^  ^ changed  from 

zero  to  one.  Otherwise  add  a Lagrange  condition  at  t to  (E*,  x,  5*?)  to  obtain 

0 p 

(Ejj,  Xjj,  ^)  • Either  way  (E^,  X^^  S?)  is  full  and  by  (2.18)  must  satisfy  the  (LPC).  By 
construction  (C)  is  satisfied  so  by  Theorem  1.1  this  is  a poised  problem.  This  is  icdpossible 
since  the  nontrivial  satisfies  the  homogeneous  problem.  Similarly  there  cannot  exist  a 

point  tj^  £ where  ♦ changes  sign  but  >(i,  does  not.  Thus  Claim  2 is 

established. 

Then 

b b 

<2-19)  /^((ij(t)i|/(t)dt  - /^|i(i*(t)  |dt  >0  , i - 1 r . 
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This  together  with  (2.17)  gives  a spline  which  contradicts  (2.13)  Thus  Claia  1 is 

established. 

Using  Leinna  2.1,  (2.14),  and  (2. IS)  we  have  that  g^  - g^  satisfies  the  )>oa»geneouB 

poised  problem  problem  (E. , X,,  S?)  and  g,  - g,  £ S*  C S?.  Thus  g,  ■ g,.  Since  g, 

Ilp  I2pp  12  1 

and  g were  arbitrary  elements  of  P_(f) , the  proof  of  the  theorem  is  complete. 

2 G 

3.  Monotonicity  and  Convexity 

Definition  (2.1)  for  G in  the  previous  section  is  a natural  generalization  to  splines 

of  the  notion  of  monotone  polynomials  introduced  by  O.  Shisha  [12]  (see  also  [4]).  Included 

is  the  possibility  for  requiring  nonnegativity  or  nonpositivity  by  choosing  • 0.  Since 

we  made  the  assumption  that  all  of  the  elements  of  our  spline  space  5*  were  continuous, 

P 

choosing  some  )(^  - 1 requires  the  usual  monotonicities,  either  nondecreasxng  or  non- 
increasing. 

If  < m - 1,  V “ l,...,q,  then  some  « 2 in^lies  that  all  of  the  elements  of 
G are  either  convex  or  concave.  However  it  is  reasonable  to  ask  for  convexity  c.  concavity 
even  if  some  of  the  knots  have  multiplicity  m - 1.  Convexity  is  well  defined  (although  not 
in  terms  of  the  second  derivative)  for  linear  splines,  i.e.  continuous  piecewise  linear 
functions,  for  example.  Similarly  monotonicity  is  well  defined  for  discontinuous  splines. 

We  briefly  indicate  how  the  preceding  section  would  need  to  be  modified  to  include 
the  requirement  of  convexity  when  ” m - 1 for  some  of  v “ l,...,q.  We  ask  that 

(2) 

(3.1)  i(i  (t)  0 , a £ t jc  b,  and 

(3.2)  ^0  whenever  « m - 1 

The  conditions  in  (3.2)  are  also  linear  constraints  on  5"*-  It  is  not  difficult  to  sliow 

P 

that  there  exists  a spline  in  which  satisfies  all  of  the  constraints  including  these 

P 

strictly  so  that  the  theorem  of  Rozema  and  Smith  110,  Theorem  4.1)  applies. 

If  none  of  the  constraints  of  type  (3.2)  are  chosen  by  the  theorem  of  Itozema  and 
Smith,  then  we  proceed  exactly  as  in  the  previous  section.  On  the  other  hand  if  one  of 
these  "jump"  constraints  is  active  and  is  chosen,  that  in^lies  that  for  g £ 
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l.e.  the  knot  \ is  really  only  of  multiplicity  m - 2 for  all  splines  in  P_(f).  It  is 

G 

easy  to  show  that  P_(f)  C p (f)  where  S®.  is  5®  with  the  knots  chosen  in  (3.3) 

■ ^ cns®  P P 

P* 

I having  multiplicity  only  m - 2 so  that  p*  < p.  In  fact  the  above  inclusion  is  an  equality 

^ and  P (f)  can  be  characterized  using  the  arguments  of  the  previous  section.  In  parti- 

' P 

t 

i cular  we  can  still  conclude  that  uniqueness  holds. 

I 

14 . Further  Extensions 

With  only  minor  modifications  the  work  of  this  paper  can  be  extended  to  the  problem 
I of  finding  a best  global  L^-approximation  to  a compact  (in  L^(a,b))  set  of  continuous 

function  F from  G as  defined  in  (2.1).  Such  best  global  approximations  are  also  called 
restricted  Chebyshev  centers  for  F with  respect  to  G or  best  approximations  to  the 
elements  of  F simultaneously.  The  methods  of  the  paper  by  Rozema  and  Smith  (10)  apply  in 
a straight-forward  manner. 

The  techniques  we  have  used  can  also  be  applied  to  the  more  general  problem  of  best 
Lj^-approximation  by  splines  with  restricted  ranges  of  their  derivatives.  In  the  uniform 
norm,  the  corresponding  polynomial  problem  was  introduced  by  Roulier  and  Taylor  [9]  and 
the  spline  problem  was  studied  in  [6] . We  wish  to  point  out  the  significant  differences 
between  uniform  and  Lj^-approximation  by  these  restricted  splines. 

Examining  the  proofs  in  [9)  and  [6J , only  the  functions  which  bound  the  ranges  of 
the  derivatives  (other  than  the  zero  derivative)  need  to  Ije  assumed  to  be  differentiable 
in  order  to  guarantee  uniqueness  in  the  uniform  norm.  In  Lj^,  the  functions  which  bound 
the  range  need  to  be  differentiable  as  well  in  order  to  carry  out  the  analog  of  Claim  2 in 
the  proof  of  Theorem  2.1.  The  problem  of  one-sided  Lj^-approximation  of  a differentiable 
function  by  splines  which  was  studied  by  A.  Pinkus  [8]  is  a special  case  where  the  given 
function  is  also  the  range  (sound. 

If  there  are  Isounds  on  the  range,  then  for  uniqueness  in  the  uniform  norm  to  be 
assured,  the  assumption  is  needed  that  the  given  function  satisfies  these  range  bounds  at 


-15- 


least  within  some  c > 0,  where  c is  strictly  less  than  the  distance  from  the  qiven 
function  to  the  set  of  restricted  splines.  If  this  is  not  the  case,  then  a single  linear 
functional  (a  point  evaluation)  in  C*[a,bl  may  be  a positive  error-extremal  and  a 
negative  constraint-extremal  or  vice  versa  in  the  terminology  of  (61 . No  such  assumption 
is  needed  in  although  the  assumption  that  the  given  function  is  continuous  is  needed. 
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